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We examine the role of magnetic dipoles in single and coupled pairs of metallic split-ring res- 
onators by numerically computing their magnitude and examining their relative contributions to 
the scattering cross section. We demonstrate that magnetic dipoles can strongly influence the scat- 
tering cross section along particular directions. It is also found that the magnetic dipole parallel 
to the incident magnetic field and/or high-order multipoles may play a significant role in the linear 
response of coupled split-ring resonators. 
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In a seminal paper published in 1999, Pendry et al. proposed that an array of metallic split-ring resonators 
(SRRs) can, in the long-wavelength approximation, possess an effective negative permeability at its inductor-capacitor 
resonance [1]. This proposal lead to the development of an artificial structured medium, a metamaterial, consisting 
' of interspaced SRR and metallic rods which simultaneously displays negative permeability and permittivity Q . Such 
novel metamaterials are absent in nature and have potential application such as subwavelength or even perfect imaging 
@j S HI • Metallic SRR arrays are a popular building block for metamaterials and have attracted intensive attention 
in the last decade 0, IE] ■ Such arrays are sometimes referred to as magnetic metamaterials due to the significant role 
. ' played by magnetic dipoles in their optical response. 

Experimentally single- and multi- layered metallic single-slit SRR membranes have been fabricated and their optical 
iT^i responses are usually studied by measuring their transmission spectra at normal incidence @, S @| • A magnetic dipole 
moment perpendicular to the membrane, arising from a current circulating inside the SRR, is frequently mentioned 
as having a role in the observed transmission spectrum and in the coupling between adjacent SRRs 0, 0, [ToL [ll[. 
However we propose that the role of this magnetic dipole is unclear or even over emphasized in the literature perhaps 
due to the difficulty in experimentally measuring magnetic dipoles. Furthermore, rarely mentioned is the contribution 
from a magnetic dipole parallel to the incident magnetic field arising from the strongly localized conduction electrons 
O I near the metallic surfaces [HJ. 

The main aim of this article is to clarify the role of these magnetic dipoles by numerically computing their magnitude 
and examining their relative contributions to the scattering cross section. To achieve our purposes, we simulate the 
optical response of an isolated single SRR and a coupled pair of SRRs using the finite-difference time-domain (FDTD) 
method Our numerical results show that the strength of the magnetic dipoles can be determined by measuring 
t-H , the scattering cross section along particular directions. We also note that in certain circumstances the magnetic 
CnJ ■ dipole parallel to the membrane plays a significant role. For isolated pairs of coupled SRRs high-order multipoles 
contribute considerately to the scattering cross section. Finally we numerically verify the qualitative assumption in 
earlier literature Q that the magnetic dipole-dipole interaction must be considered to interpret the coupling strength 
t— I ' between coupled SRRs. 

0^ . The papers is arranged as follows. In Section II we show that in the presence of an electric dipole the analytical 
expression for a magnetic dipole is not translational invariant. In section III we numerically compute the various 
electric and magnetic dipole moments and present a quantitative analysis of their relative contributions to the scat- 
tering cross section. The role of magnetic dipoles and higher order multipoles in second harmonic generation from an 
individual SRR are studied in an appendix. 



II. MULTIPOLE THEORY IN ELECTROMAGNETIC SCATTERING AND RADIATION 



Let us first recall basic multipole theory in electromagnetic scattering and radiation [12j. For a localized current 
source J, with sinusodial time dependence, the relevant vector potential at position r is given by 




where k = uj/c is the wave number, and the integration is performed over the current source. Its asymptotic form 

lim A(r) Pa J^ e lkr f J(r')e~ ifcnr V = -^_e lfcr p (2) 
j — >oo Airr J v 4irr 
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can be achieved by approximating |r — r'| as r in the denominator while in the numerator |r — r'| is replaced by r — nr' 
with n = r/r being the observation direction. Consequently, the far-zone electromagnetic fields can be formulated as 

= ik R v elkrn x E = ^ x n 
47rr 

We can now write the time-averaged power scattered per unit solid angle as 

g = I R e [r *n.E.xH;]=0|(nxp)x„|* <"» 



with rjQ = y/^to/eo being the intrinsic impedance of free space. In addition, according to the optical theorem the 
extinction power P ex taken from an incident wave Eoe* k °' r takes the form [l5[ 



P ex = ^Re 



E* • f J(r')e- lko r 'dr' 



Re [E* • p(ko)] . (5) 



This expression sug gests that P ex can be interpreted as interference between the incident wave and the forward 
scattering wave [lfj|. 

The derivations above clearly demonstrate that the direction-dependent vector p contains all the information for 
scattering and/or radiation. Using Taylor expansion, we can reformulate p as 



P = E^r : /j(r')(n-r'fdr'. 

71 ■ Jv 
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Mathematically this transformation is always valid, however, to achieve physically meaningful quantities the spatial 
dimensions of the current volume must be much smaller than the incident wavelength. In other words, the quantity 
k (n • r') should be small compared to unity, and consequently the successive terms in the expansion fall off rapidly 
with n. The first term represents the electric dipole moment, 

d = / j(r')dr', (7) 

Jv 

which corresponds to the electrostatic limit. The antisymmetric part of the second term is given by 

ikn x—fr'x j(r')dr' = ikn x m, (8) 
2 J v 

with m standing for the magnetic dipole moment, while the symmetric counterpart corresponds to electric quadrupole 

m. 

Despite the fact that these definitions can be found from any popular undergraduate electromagnetism textbook, 
determining the electric and magnetic dipole of a realistic system is not straightforward. The main reason is that as 
long as d exists m is not translational invariant (independent of the choice of origin of coordinates). More specifically 
rewriting the above in a new coordinate system x, which connects to the original as x = r + a, we obtain an identical 
far-zone vector potential A but a slightly different p given by 



p = f J(r')e- lfcn (r ' +a) dr'. 

J V 



(9) 



In this new coordinate system the electric dipole d is identical to that defined in Equation however, the magnetic 
dipole m must be reformulated as 



- / r' x J(r')dr' + x d. (10) 
2 Jv A 



The new formulation of the magnetic dipole m depends on a and is therefore not translationaly invariant, similarly 
for the high-order multipoles [13, EH • In order to compare the relative contributions of the magnetic dipoles to the 
linear scattering we therefore take the advantage of the structural symmetries by choosing the center of mass as the 
coordinate origin. 

We apply this analysis to the configurations of gold SRRs embedded in vacuum as shown in Figure 1. The structures 
are excited by a y-polarized plane wave propagating along the x direction. It excites not only a y-component electric 
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FIG. 1: Schematic of (a) an individual split-ring resonator and two coupled split-ring resonators arranged in (b) Parallel 
configuration and (c) Opposite configuration. The separation between si and s2 is 30 nm. All other dimensions shown are in 
nanometers. 

dipole d y , but also a magnetic dipole possessing both z and x components. The z component m z is due to the strongly 
localized conduction electrons of the metal, while m x is induced by the structural symmetry breaking. In order to 
simplify our analysis we ignore higher order multipoles and only consider contributions from the three dipoles d y , m x 
and m z . 

From Equation U the scattering cross section along the three primary axes are given by 

\d y =F ikm z \ 2 , if n along the ±cc direction; 

k 2 (\m x \ 2 + |m 2 | 2 ) , if n along the ±y direction; (11) 
\d y ± ikm x \ 2 , if n along the ±z direction 

For example, along the ±cc direction the forward- or backward- scattering strength is influenced by the interference 
between the d y and m z dipoles, while the ^-component magnetic dipole m x does not contribute. This observation 
is at odds with the commonly stated assumption that the m x dipole plays a significant role in one valley of the 
transmission spectrum of SRR arrays where it is often referred to as the magnetic resonance @ . Similarly, along the 
y direction the scattering cross section is due solely to the m x and m z dipoles and receives no contribution from the 
electric dipole d y . This observation may provide an insight on how to determine the relative magnitude of the various 
magnetic dipoles by experimentally measuring the scattering cross section along the orthogonal axes. 



III. RESULTS AND DISCUSSIONS 



To quantify the relative contributions of the various electric and magnetic dipole moments we perform a numerical 
simulation of the structures shown in Figure Q] using the FDTD method. The relative dielectric constant of gold is 
fitted by the Drude model of 

Ll> 2 

e(w) = l.O- * , (12) 
ujyuj + ij) 

where u> p = 1.367 x 10 16 s _1 is the bulk plasma frequency determined by the density of conduction electrons, 7 = 
6.478 x 10 13 s _1 represents the phenomenological damping rate Using FDTD, with a center of mass coordinate 
system, we compute the polarization current inside the metallic structures and employ Equations [7] and [5] to determine 
the magnitude of the electric and magnetic dipole moments. Using Equations [4] and [5l which are independent of the 
chosen coordinate system, we compute the scattering and extinction cross sections [14l |. 

As before the structures are excited with a plane wave propagating along the positive x direction and polarized in 
the y direction. The absolute extinction cross section cr cx t of the various SRR structures are shown in Figure (2a). For 
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FIG. 2: (a) Absolute extinction cross sections <r cx t. (b) Differential scattering cross section a BC , with dimensions of /im 2 per 
unit solid angle, as a function of the azimuthal angle cj>. For the sake of clarity, the curves of 1440 nm wavelength and 1540 
nm wavelength are amplified 10 times. The solid black curve corresponds to the individual split-ring resonator at 1264 nm 
wavelength. The incident plane wave propagates along the positive x [<f> = 0] direction and is j/-polarized [<j> = 90]. 



the single SRR (solid line) a peak corresponding to the fundamental plasmonic resonance was found at a wavelength 
of 1264 nm. A similar individual SRR is experimentally studied in Ref.(2fjl] and the experimental extinction cross 
section and our numerical results are in excellent agreement. A periodic array of SRRs with lattice spacing of 305 nm 
was shown to have a fundamental resonance at 1212 nm wavelength, which is very close to the isolated SSR resonance 
indicating a weak coupling between adjacent SRRs [2lj . 

The scattering cross section a sc of the isolated SRR, at 1264 nm wavelength, as a function of the azimuthal angle 4> 
in the xy plane is shown in Figure (2b) (solid line). This dependence curve can be fitted by cos 2 <fi which suggests that 
the linear scattering of this individual SRR is dominated by a y-component electric dipole moment d y . Based on the 
analysis related to Equation [IT] we can further estimate the magnitude of the magnetic dipoles from the scattering 
along the x and y directions. The plotted values of a sc along ±x directions [9 = 0, 180] are almost identical indicting 
a negligible z-component magnetic dipole m z . Similarly for scattering along the iy directions [9 = 90, 270] the cross 
section is nonzero indicating the presence of a finite total magnetic dipole. Consequently we can deduce that the SRR 
possesses a considerable x-component magnetic dipole m x . 

To compare the relative strength of the numerically computed electric and magnetic dipole moments we normalize 
them to the electric dipole moment d s of the single SRR as follows 

a=\d y /d s \, (3 X = k\m x /d s \ , /3 Z = k\m z /d s \ , (13) 

and the results are shown in Table |TJ It is found that, for the individual SRR, (3 X = 0.25 while f3 z is almost zero, 
in agreement with the scattering-cross-section-based analysis above. For comparison we also plot in Figure (2a) the 
extinction cross section of a gold square of size 220 x 220 nm. In the wavelength range of 1000 nm to 1600 nm the 
square exhibits no resonance peaks and has a weak extinction cross section consistent with its small normalized electric 
dipole of 0.23. However, it exhibits a larger z-component magnetic dipole m z than the single SRR and, because of 
its symmetry, does not present a ^-component magnetic dipole m x . 

We proceed to coupled SRRs oriented as depicted in Figure (lb,lc) with a separation between si and s2 of 30 nm. 
Two different configurations are studied. The two SRRs are parallel to each other in the first arrangement (named 
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TABLE I: Normalized magnitudes of electric and magnetic dipoles 



Structure 


Wavelength [nm] 


a[si, s 2 ] 


f3 x [si,s 2 ] 




Single 


1264 


1 


0.25 


< 0.001 


Square 


1264 


0.23 





0.001 


Parallel 


1440 


0.45 [1.65, -0.66] 


0.08 [1.93, -0.93] 


0.08 [0.68, 0.32] 


Parallel 


1208 


1.02 [0.54, 0.46] 


0.26 [0.56, 0.44] 


0.01 [4.3, -3.3] 


Opposite 


1540 


0.31 [2.06, -1.07] 


0.18 [0.58, 0.42] 


0.07 [0.63, 0.37] 


Opposite 


1176 


0.99 [0.53, 0.47] 


0.02 [9.1, -8.1] 


0.008 [4.9, -3.9] 



Parallel) while in the second arrangement s2 is rotated 180 degree in the yz plane (named Opposite). Their extinction 
spectra are shown in Figure (2a). For both configurations the 1264 nm peak of the individual SRR is split into two 
peaks induced by inter-SRR coupling. For the Parallel structure the peaks are located at 1208 nm and 1440 nm, 
and for the Opposite structure they are lactated at 1176 nm and 1540 nm. The larger separation between the peaks 
of the Opposite configuration indicates stronger inter-SRR coupling (22j. Furthermore, for each structure the longer 
wavelength peak has a magnitude comparable to that of the individual SRR, while the shorter wavelength peak has 
significantly reduced extinction cross sections. 

For each configuration the extinction cross section as a function of the angle <p in the xy plane is shown in Figure (2b). 
The scattering cross section along the y direction indicates the strength of the total magnetic dipole moment. Along 
the y direction [<f> = 90, 270] the Parallel configuration at 1208 nm possesses a scattering cross section comparable 
to that of the single SRR. This is confirmed in Table Q] where fix is 0.25 and 0.26 for the single SRR and Parallel 
configuration, respectively. For the Opposite configuration at 1176 nm, however, scattering along the y direction is 
almost zero indicating a negligible total magnetic dipole. This is also confirmed in the numerically computed values 
where (3 Z is 0.008 and (3 X is an order of magnitude smaller than the single SRR with a value of 0.02. 

The role of the z component of the magnetic dipole can be determined by looking at the difference in scattering 
between the positive and negative x directions. For example the Opposite configuration at 1540 nm has a relative 
scattering difference of 14% along these two directions. This difference is induced by interference between a y- 
component electric dipole d y and a z-component magnetic dipole m z . Consequently, we can infer that this structure 
possesses a significant z-component magnetic dipole relative to the electric dipole. Again this is confirmed in Table [I] 
where the ratio of magnetic to electric dipole strength j3 z /a has a value of 0.23. 

Higher order multipoles play a more profound role in the linear scattering of coupled SRRs since their effective 
thickness of 80 nm is much larger than 25 nm of a single SRR. At the shorter wavelength peaks both coupled 
configurations possess an electric dipole comparable to that of the single SRR, however, their scattering along the 
x direction [9 = 0] differ by 12%. Taking into account that (3 Z is negligible for all three structures this difference 
cannot be interpreted by interference between the electric and magnetic dipoles alone. Consequently, contributions 
from higher-order multipoles, such as electric quadrupole, must be included. It should be mentioned that higher-order 
multipoles are also found to affect the optical properties of stereometamaterials Q- 

To distinguish the relative contribution from each SRR in the coupled structures we compute their individual dipoles 
and normalize them to the total dipole moment. These values are listed in square brackets in Table [J where values 
with opposite sign indicate the relevant dipoles are out of phase. For both Opposite and Parallel configurations the 
y-component electric dipoles d y are in phase at the shorter wavelength peaks and out of phase at longer wavelength 
peaks. This is consistent with the fact that a coupled system possesses a high-energy symmetric mode as well as 
a lower-energy anti-symmetric mode. The inverse of this rule applies to the z-component magnetic dipoles m z in 
that they are out of phase at shorter wavelengths and in phase at longer wavelengths due to the different symmetry 
of electric and magnetic fields. Because of the structural symmetry a more complicated phase relationship exists 
between the x-component magnetic dipoles m x . For example, at the shorter- wavelength resonance they are in phase 
for the Parallel configuration while out of phase for the Opposite one. Table Q] further implies that although electric 
dipole-dipole interaction dominates the inter-SRR coupling, magnetic dipole-dipole interaction and in particular the 
m x — m x interaction must be included to interpret the fact that the Opposite configuration presents a larger separation 
between resonance peaks than its counterpart. 

IV. CONCLUSIONS 

This paper clarifies the role of magnetic dipoles in single and double metallic split-ring resonators (SRRs) by 
numerically computing their magnitudes and examining their relative contributions to the scattering cross section. 
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FIG. 3: Linear (a) and second-harmonic (b) scattering of the individual split-ring resonator, specified in Figure (la), in the yz 
plane. Both y-polarized (90 degree) and z-polarized (zero degree) incidences are considered. As a reference, a metallic square 
is also studied. 

We propose that the strength of the magnetic dipoles can be determined by measuring the scattering cross section 
along particular directions. It is demonstrated that the magnetic dipole parallel to the incident magnetic field plays 
a significant role for Opposite configuration at its longer wavelength resonance. Moreover, we found that for isolated 
pairs of coupled SRRs higher-order multipoles contribute considerately to the scattering cross section. By considering 
the relative phase between the partial dipoles possessed by individual SRRs, in coupled configurations, we have shown 
that the electric dipole-dipole interaction dominates the coupling between adjacent SRRs and also that the magnetic 
dipole-dipole interaction is responsible for the larger resonance separation in the Opposite configuration. 

This work is supported by the Air Force Office of Scientific Research (AFOSR), under Grant No. FA9550-07- 1-0010 
and FA9550-04- 1-0213. J. V. Moloney acknowledges support from the Alexander von Humboldt foundation. 



Using a recently developed classical theory [231 ], we numerically studied second-harmonic radiation from an indi- 
vidual SRR as schematically shown in Figure (la). The fundamental plane wave propagates along the positive x 
direction and has a wavelength of 1264 nm, corresponding to the fundamental resonances excited by a y-polarized 
light. Its linear scattering cross section in the yz plane is shown in Figure (3a). For comparison we also considered a 
gold square of size 220 x 220 nm. Its linear response is quite close to that of the SRR with a z-polarized excitation, 
that is, they both possess an almost identical-strength electric dipole and a negligible magnetic dipole. The strongest 
electric dipole appears in the SRR with a y-polarized excitation. 

Figure (3b) plots second-harmonic radiation in the yz plane. Because of structural symmetry the gold square 
presents a dominant electric quadrupole pattern. On the other hand, magnetic dipoles are found to contribute 
considerably to the SRR with the y-polarized excitation, and the radiation exhibits a slightly distorted dipole pattern. 
It is mainly due to the fact that the second-harmonic wavelength is around 600 nm while the lateral size of the SRR 
is roughly 200 nm. In addition, although the square has almost identical linear response as the SRR with z-polarized 
excitation, its second-harmonic strength is much weaker. This is because its radiation is dominated by an electric 
quadrupole, while the SRR radiation is dominated by the electric and magnetic dipole. 
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